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We demonstrate a fast production of large 23 Na Bose-Einstein condensates in an optically plugged, 
magnetic quadrupole trap. A single global minimum of the trapping potential is generated by slightly 
displacing the plug beam from the center of the quadrupole field. With a dark magneto-optical trap 
and a simple rf evaporation, our system produces a condensate with N « 10 7 atoms every 17 s. 
The Majorana loss rates and the resultant heating rates for various temperatures are measured with 
and without plugging. The average energy of a spin- flipped atom is almost linearly proportional to 
temperature and determined to be about 60% of the average energy of a trapped atom. We present 
a numerical study of the evaporation dynamics in a plugged linear trap. 

PACS numbers: 03.75.Hh, 37.10.De 



I. INTRODUCTION 

Ultracold quantum gases have presented a new class 
of systems to investigate many-particle quantum phe- 
nomena [TH3]. Since ultracold gas samples are very di- 
lute at particle densities around 10 13 or 10 14 /cm 3 , we 
need extremely low nanokelvin temperatures to reach 
quantum degenerate regime. Evaporative cooling has 
been the powerful and exclusive method for producing 
quantum gas samples such as Bose-Einstein condensates 
(BECs) [4]-[6]. Rapid and efficient evaporation requires 
fast rethermalization and conservative environment min- 
imizing undesired atom loss and heating due to, e.g., in- 
elastic collisions and background collisions [7 J. In this 
paper, we describe a fast evaporation to produce a large 
condensate in an optically plugged magnetic quadrupole 
trap. 

A quadrupole trap provides a large trap volume and 
tight confinement due to the linearity of the potential, 
which make it beneficial to efficient evaporation. Also 
the quadruple field cab be generated from a simple pair 
of coaxial coils with opposite currents, allowing large op- 
tical access to a sample. However, Majorana spin- flip loss 
is a significant drawback in the quadrupole trap. Near 
the zero- field point at the trap center, nonadiabatic spin 
transitions to an untrappable state cause atom loss and 
resultant heating, eventually limiting the maximum gain 
in phase-space density [8 j. Several techniques have been 
implemented to avoid this detrimental effect, including 
applying a time-orbiting bias field [4j or an repulsive op- 
tical plug potential [5] |9], or transferring to an optical 
dipole trap [TO] . 

The optical plugging scheme preserves the linearity 
of the trapping potential, so it is advantageous for fast 
evaporation. The first 23 Na BEC of « 5 x 10 5 atoms 
was produced in 7 s evaporation [5 j. Recently, high 
efficiency of evaporation cooling in a plugged trap has 
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been demonstrated by producing a large condensate of 
w 3.0 x 10 T atoms in 40 s evaporation [9 j. Here we 
report a rapid production of a large condensate with 
N « 1.0 x 10 T atoms in 14 s evaporation. The distinct 
feature of our system is that the plug beam is offsetted 
from the magnetic zero-field point. This seems contra- 
dicting to the role of the plug, but we find that a slight 
misalignment less than a beam waist is favorable for BEC 
production because it forms a single global minimum in 
the combined magnetic and optical potential. The point- 
ing stability of the plug beam was improved by using a 
fiber, which makes the fine adjustment of the plug posi- 
tion possible. 

Furthermore, we have investigated the effects of the 
Majorana loss and the plug by measuring the loss rates 
and the resultant heating rates for various conditions. 
The average loss energy of a spin-flipped atom shows an 
almost linear dependence on temperature and is deter- 
mined to be about 60% of the average energy of a trapped 
atom. For a plugged trap, the reduction factor in the 
loss rate is found to be exponentially proportional to a 
Boltzmann factor. This confirms that the plugging ef- 
fect results from density suppression due to the repulsive 
optical potential in the loss region. Finally, we present 
numerical simulation for quantitative description of the 
evaporation dynamics in a plugged quadrupole trap. 



II. EXPERIMENT 

Our BEC apparatus employs a zero crossing Zee- 
man slower, designed to generate a flux of about 10 11 
atoms/s [TT] , Slowed atoms were collected for 3 s in a 
dark magneto-optical trap (MOT) [12] and the atoms 
in the hyperfine state \F = = —1) were trans- 

ferred in a magnetic quadrupole trap, where F is the 
total angular momentum and mp is the magnetic quan- 
tum number. The magnetic field was generated from a 
pair of 36 turn coils, secured around recessed bucket win- 
dows of our main experimental vacuum chamber [Fig. [I] . 
The coils are made from 1/8 in. hollow square copper 
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FIG. 1: Experimental apparatus. Recessed bucket windows 
with 60 mm diameter clearance are mounted above and below 
the main vacuum chamber. A pair of 36 turn water-cooled 
coils are secured around the windows, generating a magnetic 
quadrupole field at the chamber center. Radio-frequency (rf) 
evaporation coils are placed inside the vacuum chamber. A 
plug beam is aligned with the axial direction of the magnetic 
field. 



tubes and cooled with water flowing through the tubes. 
The average coil diameter is 90 mm and their spacing is 
about 63 mm. We measured a gradient of 1.5 G/cm/A 
along the symmetry z-axis. For the MOT and the mag- 
netic trap, we applied a field gradient B' = 7.5 G/cm 
and 180 G/cm with 5 A and 120 A in the trap coils, 
respectively. The atom transferring was simply done by 
ramping B' at a rate of about 10 G/cm/ms and simul- 
taneously turning off the MOT laser beams. The initial 
number and temperature of the atoms in the magnetic 
trap were Ni ~ 3.0 x 10 9 atoms and Ti ~ 250 /iK, re- 
spectively. 

Optical plugging was performed by focusing a blue- 
detuned, 532 nm laser beam near the zero of the magnetic 
field. The plug beam propagated along the z-direction 
with a beam waist of w ~ 45 /im, providing a repulsive 
potential barrier of U p c± ks x 90 /iK with 1.7 W, where 
fc# is the Boltzmann constant. In order to improve its 
pointing stability, we used a fiber to deliver the plug laser 
beam close to the main chamber, which was critically im- 
portant for our BEC production [13]. For several weeks 
of operation, the plug position needed little or no adjust- 
ment. Finally, rf evaporative cooling was applied for 14 s 
to the atoms in the plugged quadrupole trap. 



A. Quadrupole trap with off-centered plug beam 

When a circular plug beam is placed at the center 
of the quadrupole field, the combined potential has az- 
imuthal symmetry, forming a ring trap [Figj2fa)]. This 
ring geometry is interesting in terms of persistent cur- 
rent and vortex generation [14]. However, it is desirable 
to break this symmetry and generate a single point min- 
imum in the trapping potential for controlled production 
of condensates. Furthermore, the critical temperature in 
a ring trap for a given atom number is typically low be- 
cause of no confinement along the azimuthal direction. 
In previous experiments, this ring symmetry was broken 
by applying a plug beam in the transverse direction [5] or 
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FIG. 2: (Color online) Off-centered plugged trap. Schematic 
diagrams of the combined trapping potential with a plug 
beam (a) passing the field center and (b) slightly displaced 
from the center. The ring symmetry is broken by the displace- 
ment and a single minimum is generated in the trap. In-situ 
absorption images of atoms in a ring trap after evaporation 
with a final rf of (c) 150 kHz and (d) 110 kHz, where multiple 
small condensates containing < 10 6 atoms are produced due 
to the imperfection of the plug beam profile, (e) In-situ image 
with displacing the plug beam by d ~ 17 /im and a final rf 
of 105 kHz. A condensate of ~ 10 7 atoms is produced in a 
single potential minimum. The white trace lines in the im- 
ages are from the leakage plug beam out of filters, which are 
used for determining the relative beam position. The field of 
view is 750 /mix 750 /im. (f) Radial potential along the trap 
minimum at B' = 30 G/cm (solid) and 180 G/cm (dashed). 



using an elliptical beam [9 . In our experiment, we simply 
displaced the plug beam from the zero of the magnetic 
field by ~ 17 /im [Fig. [2^b)] . It has been emphasized 
that misalignment of a plug is detrimental to BEC pro- 
duction [9] . We found that a slight misalignment smaller 
than a beam waist is effectively helpful without signifi- 
cantly losing the plugging effect. The position stability 
of our plug beam was estimated to be less than 5 /im. 
The effective potential for atoms in the plugged trap 

is 



U(r) = iiB' 



4 



(x + df 
Up exp( ^2—) - m 9 z ^ 

(1) 

where \i = /i#/2 (/i# is the Bohr magneton) and m are 
the magnetic moment and mass of the atom, respectively, 
d is the plug beam offset from the field center, and g is the 
gravity acceleration. We have ignored the beam focusing 
because the Rayleigh length ttw 2 /\ « 12 mm is larger 
than the sample size. For our typical final condition with 
B' = 30 G/cm and d = 17 /im, we estimate the trap min- 
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FIG. 3: Experimental sequence for BEC production consists 
of dark MOT, rf evaporation and trap decompression, (a) 
Field gradient and radio-frequency, (b) Time evolution of 
atom number, temperature, truncation parameter 77, peak 
density, and collision time. The decompression phase is 
marked with a grey zone. 



imum to be located at (xo,^o,^o) = (60,0,15) jam with 
a magnetic field of about 100 mG. From the potential lo- 
cal curvature, the trapping frequencies are estimated to 
be (fxjyjz) ~ (210,55,110) Hz. With d = OAw, the 
plugging potential at the zero-field point is reduced by 
about 25% to ks x 68 /iK. At low temperatures, we ob- 
served multiple atom clouds due to small imperfections 
in the profile of our plug beam [Fig.^d)]. The potential 
variation leading to this fragmentation was determined 
to be about hxlO kHz from the density distributions for 
various final rf values, where h is the Planck constant. 
In the off-centered plugged trap, a single potential mini- 
mum was formed at h x 30 kHz lower than next potential 
minimum [Fig. [2^e)]. 



reducing inelastic losses associated with high atom den- 
sity such as three-body recombinations. Furthermore, 
decompression results in adiabatic cooling of atom cloud, 
thus relatively enhancing the plugging effect. During the 
adiabatic process in a linear trap, atom density and tem- 
perature scale as B' and £> /2 / 3 , respectively. 

The time evolution of atom number and temperature 
during the evaporation was determined from absorption 
images taken after 18 ms time-of- flight [Fig. 3j. Peak 
density n and collision time r are estimated as n = 
N/(327r(k B T/fiB') 3 ) and r -1 = V2nav, respectively, ig- 
noring gravity and the plug potential, a = 87m 2 is the 
elastic collision cross section with scattering length a — 
52an (an is the Bohr radius) [15] and v = (SA^T/Vm) 1 / 2 
is the average thermal velocity. The truncation param- 
eter r] is defined as 77 = (hv x { — Uo)/kBT, where Uq is 
the potential at the trap bottom, Uo/h ~ 320 kHz at 
180 G/cm and « 70 kHz at 30 G/cm [Fig. ^f)]. Run- 
away evaporation was achieved with increasing collision 
rate r _1 . When the decompression started, the peak 
density was estimated to be ~ 10 14 /cm 3 . rj smoothly 
increased from 6 to 9 before the decompression. 

Bose-Einstein condensation was observed after 13.7 s 
evaporation from the onset of bimodal distribution in ab- 
sorption image [Fig. Qb)]. At this transition point, the 
atom number and the temperature were N c = 2.5 x 10 7 
and T c = 1.4 /iK, respectively. The overall evapora- 
tion efficiency is 7 = -\n(D c /Di)/ \n(N c /Ni) = 2.7, 
where D is peak phase space density given as D = n\^ dB 
(^dB = h/\/27rmkBT is the thermal de Broglie wave- 
length) and the critical phase space density is D c = 2.61. 
With 14 s full evaporation, an almost pure condensate 
with N = 1.0 x 10 7 was obtained. After an additional 
decompression to B' = 15 G/cm, the lifetime of the con- 
densate was 10 s, which is long enough for further manip- 
ulation such as transferring into an optical dipole trap. 

The evaporation in a bare quadrupole trap was also 
investigated without the plug beam. Figure [4] shows the 
evaporation trajectories in the N-T and N-D planes. At 
T < 20 /iK, the plugging effect becomes noticeable from 
the deviation between the evaporation curves with and 
without plugging. The maximum phase space density 
was about 5 x 10 -2 in a unplugged quadrupole trap. 
The minimum plug beam power for making a condensate 
of > 10 6 atoms was ~ 0.4 W for the 14 s evaporation 
[Fig. [5]. In our experiment setup, we could minimize the 
evaporation time to 8 s to observe a condensate. 



B. Evaporation process 

Evaporative cooling was performed by linearly sweep- 
ing the applied radio frequency v v { from 30 MHz to 2 MHz 
in 12 s and then to 0.1 MHz in another 2 s [Fig. [3ja)]. 
The quadrupole field gradient B' was ramped down from 
180 G/cm to 30 G/cm in the last 2 s. This trap de- 
compression was necessary to achieve a large condensate, 



C. Effects of Majorana loss 

The Majorana spin flips cause heating as well as atom 
number reduction. To quantitatively understand the loss 
effects, we have measured the loss and heating rate of a 
thermal gas for various temperatures with and without 
the plug at £'=180 G/cm [Fig. [6]. The initial sample 
temperature was controlled with the final value of ra- 
dio frequency. The time evolution of atom number and 
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FIG. 4: (Color online) Evaporation trajectory, (a) Temperature versus atom number, (b) Peak phase space density versus 
atom number. Data points at 11.9 s (before decompression) are marked with 'K'. The onset of condensation is observed at 
13.7 s in a plugged trap from the bimodal momentum distribution (see inset). The black diamond indicates the transition 
point. The solid lines are from num erica l simulation with a plug (black), no plug (red), and no loss (blue). The dashed lines 
are the evolution after 12 s (see Sec. III). 
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FIG. 5: (Color online) Condensate number versus plug power 
Each point is averaged over three separate runs of the experi- 
ment. Maximum phase space density (blue dashed) and con 
densate number (red) from numerical simulation (see Sec. Ill ) 



temperature were determined with various holding times. 
The rf field was turned off after the sample preparation in 
order to exclude the evaporation effects [16 j. We observed 
that the temperature almost linearly increases with a 
small acceleration and the number decay rate decreases 
slowly. The initial loss and heating rate were determined 
from an exponential and a linear fit, respectively, to the 
data points with temperatures lower than 120 % of the 
initial temperature. 

A simple argument in Ref. [17 suggests that the 
Major ana loss rate T m in a quadrupole trap scales as 
r m ex h/m(nB' /k B T) 2 . A fit to our data in Fig. [(|a) 
with CT~ 2 + T b yields C = 139(3) /iK 2 /s and T b = 
0.013(7) s _1 , where T^ is the background loss rate (cor- 
responding to the lifetime of 75 s). This result suggests 
T m = 0.14 h/m(nB' /k B T) 2 . The proportionality con- 
stant is found to be much smaller than the values of 



0.58 [T7] [18] and 0.87 [10] in previous experiments with 
87 Rb atoms, which is beyond the simple model. Note 
that the lifetime of 23 Na atoms in a quadrupole trap is 
comparable to that of 87 Rb atoms even with their mass 
difference. 

A spin flip happens when the spin of a moving atom 
cannot adiabatically follow the field orientation, so the 
flipping probability of an atom at a given position is fully 
determined by the local field strength, its gradient, and 
the velocity distribution of the atom, i.e. temperature. 
Therefore, the plugging effect is a consequence of the 
density suppression due to the optical repulsive potential 
in the low-field, spin-flipping region. Including a Boltz- 
mann density suppression factor, the reduced loss rate 
would be given as 



r^ = r m x/exp(-f/7fc B r), 



(2) 



where U* is the characteristic potential of the plug and / 
is a factor depending on the details of the plug beam such 
as shape, position and power. We find this model nicely 
fit to our data with C'T~ 2 exp(-/7* /k B T) + T 6 , where 
U* = k B x 25(1) /iK and C = 97(6) /iK 2 /s (/ w 0.7). 
This exponential behavior clearly supports our picture of 
the plugging effect. Note that U* is consistent with the 
deviation point in the N-T evaporation trajectories and 
close to the potential at the trap bottom Uo ~ k b x 15 /iK. 
In our experiment, the collision rate was approximately 
10 3 times faster than the loss rate, validating the thermal 
equilibrium assumption. 

We characterized the heating effect of the Majorana 
loss with determination of the average loss energy of a 
spin- flipped atom £ m k B T. From the energy conservation 
relation E = —e m k B TN , where E = Nek B T is the total 
energy of the system with the average energy of a trapped 
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FIG. 6: Majorana loss in a magnetic quadrupole trap, (a) 
Loss rates and (b) heating rates for various temperature at the 
field gradient B x =180 G/cm. Solid curves are fits to the data 
with a loss model, Eq. (pi (see text for detail), (c) Average 
energy of a lost atom due to the Majorana spin- flip. The dot 
line indicates the average energy of an atom in a linear trap. 



atom eksT, we derive a simple relation between the tem- 
perature and the atom number, T/T = (e m /e — 1) TV/TV. 
When e m < e, the system is heated up and vice versa. 
In a power-law trap U oc r 3 / 5 , e = (3/2 + 5). Our mea- 
surements show that Em has a very weak temperature 
dependence, slightly increasing at low temperature with 
the average value of 2.8, which is less than £=4.5 for a 
linear trap (5=3) [Fig. [6^b)]. Remarkably, we observe 
that e m is not affected by the optical plug in the temper- 
ature range of 10 to 100 /iK. This is somewhat surprising 
because the average energy of a lost atom should increase 
in the presence of the plug due to the repulsive optical 
potential. We estimate that the effective increment in the 
average loss energy would be U* — Uo « x 10 /iK with 
respect to the trap bottom and it might be observable at 
lower temperatures. 



III. NUMERICAL SIMULATION 

In this section, we numerically study the dynamics of 
evaporation process in an optically plugged quadrupole 
trap. Including the only three loss processes: background 
collision, rf evaporation and the Majorana spin- flip, the 
rate equations for atom number N and temperatue T are 



TV 
TV 
t 
f 



= -( 



IV 



i)r e -(- 



i)r m , 



(3a) 
(3b) 



where T e is the evaporation loss rate and e e = 77 + ft is the 
average energy of an evaporated atom in fc^T. The first 
term in Eq. ([3J3) corresponds to the cooling from evap- 
oration. We assume that the background collision loss 
occurs equally to all trapped atoms, thus not affecting 
the temperature. 

The simulation for our experiment was performed 
in the following conditions: TV(0)=TV; and T(0)=T* 
at t = s, i/ rf (t)=(30-7t/3) MHz, r]=hu vi /k B T, 
r5=0.013/s and £ m =0.28. The Majorana loss rate 
is parameterized with the plug beam power P as 
r m =C(l - aP)T- 2 exp(-bP/T), where C=139 /iK 2 /s, 
a=0.176/W and 6=1.48 /iK/W. For 3D evaporation in 
a linear trap (5=3), standard evaporation theory [19] 
gives Yf^ncrve-^Pd^) - ^P(^^)}/P(lv) 2 and 
k=1 - P(f ,v)/[vP(hv) - t P (t^)]> respectively, 
where P(a,rj) is the incomplete gamma function [20] . 
However, the observed evaporation speed was almost 2 
times faster than the theory value. This is quite surpris- 
ing because the theoretical estimation accounts for the 
maximum 3D evaporation, worthy of further investiga- 
tion. We empirically set Y e = 1.95 x Tf 1 in our simula- 
tion. 

The simulated evaporation trajectories are presented 
with the experiment data in Fig. [4] The trap decompres- 
sion is not included in our simulation, so the trajectories 
deviate from the data points after 12 s. The maximum 
phase space density and the critical atom number TV C 
at the phase transition were numerically estimated as a 
function of the plug beam power [Fig. [5], which show 
a good agreement with the experiment results. We ap- 
proximate the condensate number at 25% of the critical 
atom number. The simulation with a perfect plug, i.e. 
r m =0 (blue lines in Fig. [4| suggests that the maximum 
condensate would have TV w 2 x 10 7 atoms for the given 
14 s evaporation. 



IV. CONCLUSIONS 

We have described a rapid production of 23 Na Bose- 
Einstein condensates in an optically plugged magnetic 
quadrupole trap. The off- axis optical plug generates a 
single global minimum in the combined trapping poten- 
tial and consequently enhances the BEC production. The 
Majorana loss, which is critical in a quadrupole trap, has 
been quantitatively investigated, determining the aver- 
age loss energy of a spin-flipped atom and the loss rate 
suppression of the optical plug. We have presented a 
simple model to explain the evaporation dynamics in the 
plugged quadrupole trap. This model will be helpful to 
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develop and optimize evaporation processes using a mag- 
netic quadrupole trap. 
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